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We present a theoretical study on the superconductivity of /3'-(BEDT-TTF)2lCl2 at 
Tc —14.2 K under a high hydrostatic pressure recently found, which is the highest among or- 
ganic superconductors. In the present work, we study an effective model using the fluctuation- 
exchange (FLEX) approximation based on the results of first-principles calculation. In the 
obtained phase diagram, the superconductivity with d^^y-like symmetry is realized next to the 
antiferromagnetic phase, as a result of the one-dimensional to two-dimensional crossover driven 
by the pressure. 



BEDT-TTF (ET) compounds can have many kinds of 
layered structures, each of which possesses a distinctive 
function. One of the most prominent functions is super- 
conductivity, such as K-(ET)2X.^^ It is shown that the 
temperature-pressure phase diagram of k-(ET)2X salt 
can be reproduced by the fluctuation-exchange (FLEX) 
approximation,'^"^) which is a kind of self-consistent spin- 
fluctuation theory.^' According to the FLEX approxi- 
mation, d-wave superconductivity is expected and is me- 
diated by the strong antiferromagnetic (AF) fluctuations 
due to the Coulomb interaction. 

Recently, Taniguchi et al. discovered another super- 
conductor, /3'-(ET)2lCl2, whose superconducting (SC) 
transition temperature, Tc, is 14.2 K under a high hydro- 
static pressure (P >^ 8.2 GPa). This is a new record 
among organic superconductors. At ambient pressure, 
/?'-(ET)2lCl2 shows a semiconducting temperature de- 
pendence of conductivity below room temperature, and 
indicates a magnetic transition at Tn=22 K. As pressure 
is applied, resistivity decreases gradually, and metallic 
behavior [dp/dT > 0) is observed above Tmit under 
6.5 GPa. Note that Tmit, which is the crossover tem- 
perature, is expected to be higher than Tn- At 8.2 GPa, 
the insulating phase (or Tmit) disappears and the SC 
transition occurs at Tc = 14.2 K at the same time. 

At present, the atomic positions of /3'-(ET)2lCl2 under 
applied pressure are unknown experimentally. In order 
to explain the electronic property in /3'-(ET)2lCl2 under 
pressure, Kontani discussed the phase diagram of the ti- 
tle material in connection with the /3-(ET)2X, which is a 
superconductor (Tc = 1 ^ 8 K) with a two-dimensional 
Fermi surface. He proposed that d^2_y2-wave super- 
conductivity due to AF fluctuations occurs, driven by 
the dimensional crossover under applied pressure. 

After Kontani's study, the optimized atomic positions 
of /3'-(ET)2lCl2 under appHed pressure were calculated 
using the first-principles method based on the gener- 
alized gradient approximation.^' At the same time, 
tight-binding parameters were fitted to reproduce the 
electronic structure near the Fermi level on the b*-c* 
plane at several pressures. According to the result, the 
Fermi surface possesses two-dimensional nature at higher 



pressures as expected in ref. 9. However, its shape is dif- 
ferent from that of /3'-(ET)2lCl2, contrary to the assump- 
tion in ref. 9. 

In the present work, we study the origin of the su- 
perconductivity as well as the phase diagram in /?'- 
(ET)2lCl2 based on the newly derived electronic struc- 
ture from the first-principles calculation. We employ 
the FLEX approximation based on the fact that the SC 
phase appears upon degradation of the AF phase by pres- 
sure. Although the obtained AF region is rather wider 
than that in experiments, the derived Tc is consistent 
with experiments on /3'-(ET)2lCl2. The present analysis 
suggests that dxy-wave superconductivity is realized due 
to spin fluctuations, Q = {qy,qz) ~ (tTjO). 

Although the structures of organic materials are com- 
plex, it is well known that the band structure near the 
Fermi level is quite simple; it is sufficient to take into 
account the HOMOs of the ET molecules for the title 
material. The ET network is two-dimensional, refiecting 
the layered structure, and the schematic structures of (3'- 
(ET)2X are shown in Fig. 1(a). Each ellipse represents 
an ET molecule, and each highest occupied molecular or- 
bital (HOMO) of the ET molecule possesses 1.5 electrons 
on average. The fitted transfer integrals are shown in Ta- 
ble I and are plotted in Fig. 1(b). As a good approx- 
imation, we take into account only the antibonding orbit 
of each pair of ET molecules connected by t{pl), whose 
absolute value is much larger than the others. Then, the 
original system is well mapped onto the "dimer model" 
at half filling. Hereafter, we study the dimer Hubbard 
model, whose relevant parameters are the tight binding 
integrals obtained by the fit to the upper HOMO band 
of the ET molecules, and the effective on-site Coulomb 
interaction [Ues) on the dimer molecule. 

First, let us review the electronic structure at ?7 = 0. 
In this study the pressure effect on the electronic states is 
only taken into account through the change in the trans- 
fer integrals: they are interpolated by the cubic spline 
up to a pressure of 12 GPa and extrapolated linearly 
over a pressure of 12 GPa. (See also Fig. 1(b).) When 
one compares each transfer integral at ambient pressure 
with those at 4 GPa, the transfer integrals change greatly. 



2 J. Phys. Soc. Jpn. 



Letter 



Hiori Kino, Hiroshi Kontani, and Tsuyoshi Miyazaki 




> 



/ I -0.4 





(d) 


1 ''=" — 


P=4 - - 




1 '■=8 -- 






\\\:'- 


\J.\ 1 P=16 — 















-0.3 -0.2 -0.1 0.1 0.2 Q.i 

(eV) 



(e) 






















F=0 
























P=12 









Fig. 1. (a) Schematic crystal structure of 6c-plane and transfer in- 
tegrals between BEDT-TTF molecules, (b) Fitted transfer inte- 
grals (c) band width, (d) the density of states and (e) Fermi sur- 
faces at applied pressures. The first-principles electronic struc- 
tures are calculated structures at 0, 4, 8, and 12 Gpa. The trans- 
fer integrals at other pressures are interpolated using cubic spline 
or extrapolated linearly. The unit of P (pressure) is GPa. 



However after 6 GPa, they hardly change except for the 
hnear increase in t(jj2), which is the counterpart of the 
intradimer transfer integral, t(pl). As a result, the width 
of the upper band increases almost linearly after 4 GPa. 
(Fig. 1(c)) Such a feature is also confirmed in the density 
of states shown in Fig. 1 (d) . The Fermi surfaces are also 
shown in Fig. 1(e) at a number of pressures. At an am- 
bient pressure, the Fermi surface looks one-dimensional 
and the nesting is strong. However, at higher pressures, 
the Fermi surfaces become two-dimensional and a closed 
Fermi surface appears at P=16.0 GPa, because of the 
upward shift of the saddle point near the F point in the 
electronic structure. Note that the saddle point is located 
at the r point, which results in the van Hove singular- 
ity in the density of states. It is located just below the 
Fermi level at P=12.0 GPa and above the Fermi level 
at P=16.0 GPa. (also see Fig 1(d).) This peak gives 
a larger density of states than that expected from the 
bandwidth. This point will be discussed later. 

Next, in order to take into account the effects of elec- 
tron correlations, Ucs, we calculate the self-energy for 
the dimer model for C/off=0.25, 0.30 and 0.35 eV using 
the FLEX approximation.^^ Hereafter, we write UcS as U 
for simplicity. To obtain the magnetic transition temper- 
ature Tn, we calculate the Stoner factor without vertex 
corrections, as, given by 



as = max{ ?7 • x°(k, c^ = 0) }, 



(1) 



Table I. Fitted transfer integrals. Units of pressure and transfer 
integrals are GPa and eV, respectively. ^''^ See also Fig. 1(a) to 
specify the location of the transfer integrals in the crystal. 
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Fig. 2. Neel temperature (Tn) and superconducting transition 
temperature (Tc) as functions of pressure. is evaluated by 
the Stoner criterion, as=0.995 or 0.997. Lines are guides for the 
eyes. The unit of U is eV. 



where 

X°(q, u:i) = -T^ G(q + k, c., + e„)G(k, e„), (2) 

k.n 

where G(q + k,a;/ + e„) is the Green function given by 
the FLEX approximation, and loi (e„) is the Matsub- 
ara frequency for boson (fermion). Tn is determined by 
the Stoner criterion, as = 1. In the FLEX approxima- 
tion, however, as does not exceed 1 at finite temper- 
atures in two-dimensional systems, which is consistent 
with the Mermin and Wagner theorem. So we deter- 
mine Tn under the condition as = aN, where we set 
aN as (1 — aN)~^ ^ O(IOO). The AF state will occur 
through the weak coupling between layers, Jj_.^^^ To 
evaluate the SC transition temperature Tc, We solve the 
linearized Eliashberg equation. For a singlet-pairing case 
[0(-k, e„) = -f(/)(k, e„)], it is given by^^^ 



A • (/)(k, e,, 



-T^y(k-q,e, 



xG(q, e„JG(-q, -e. 



0(q, e„)(3) 



where nk,c.O - W ' tW^) ' l''^ ' T^uW!^) + 
Tc is given by the condition A = 1. In the FLEX approx- 
imation, a finite Tc is obtained even in two-dimensional 
systems irrespective of the Hohenberg theorem. However, 
this approximation gives reasonable Tc's for k-(ET) or- 
ganic compounds and high-Tc cuprates.'^"^' In these 
calculations, we adopted 64x64 fc-points and 512 Mat- 
subara frequencies, and executed summations on fre- 
quency and wave vectors by employing fast Fourier trans- 
form techniques. 
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The obtained phase diagram is given in Fig. 2. First 
let us see the Neel temperature. At U =0.30 eV, the Tn 
at ambient pressure is 16 K. With increasing pressure, 
Tn increases rapidly and has a peak at 75 K at 6 GPa. 
With further increase in pressure, Tn decreases gradually 
and finally Tn becomes smaller than Tc at 15.2 GPa and 
the antiferromagnetism is substituted for the supercon- 
ductivity in the phase diagram. A similar trend is also 
observed at other f/s; e.g., Tn is 12 K at ambient pres- 
sure at U=0.25 eV, has a maximum value at 69 K at 
6 GPa and decreases under further applied pressure and 
is finally replaced with Tc. A larger U yields a larger Tn 
as well as a wider area of the antiferromagnetic phase. 

Next we discuss the obtained magnetic order. The 
nesting vector Q where x((5,w) has a maximum is 
commensurate and (tt, 0) for the low-pressure region, 
P <12 GPa for C/=0.35 eV (P <8 GPa for t/=0.30 and 
U =0.25 eV). This is a result of the strong electron corre- 
lation given by the FLEX approximation; in fact, Q takes 
an incommensurate value at t/ = 0. For higher- pressure 
regions, on the other hand, Q becomes incommensurate 
and is (tt — 0.1257r, 0) when Tc has a maximum. A shallow 
dip in Tn at approximately 14 GPa for C/=0.35 (10 GPa 
for J7=0.25 and [/=0.30) corresponds to the crossover 
pressure between the commensurate and incommensu- 
rate Q. 

Finally let us discuss the superconductivity. If the 
Stoner criterion is 0.995, the maximum Tc under the 
condition Tc > Tn for U = 0.30 eV is approximately 
18 K at 15.2 GPa. It decreases as pressure increases, 
and Tc ~ 8 K at 17.6 GPa. The same feature can be 
seen in the case of C/=0.25 eV and 0.35 eV except for 
the shift in the transition temperatures. A smaller U 
gives a smaller Tc as well as a smaller transition pressure 
between antiferromagnetism and superconductivity. The 
maximum Tc is 9.5 K(18.5 K) at 14.4 GPa(17.6 GPa) for 
C/=0.25 eV(0.35 eV). 

Figure 3 shows the solution of eq. (3) at (fcy, kz) on the 
Fermi surface. The inset shows the change of the Fermi 
surface at C/=0 and [7=0.30 eV. The saddle point near 
the F point causes a peak in the DOS above the Fermi 
level at C/=0, however, it is below the Fermi level at 
U =0.3 eV because of the fc-dependence of the self-energy 
(Sfc(O)). As a result, the Fermi surface for t/=0.3 eV is 
still open at this pressure. 

The main figure shows that the obtained SC order pa- 
rameter has line nodes (along a direction) as in the high- 
Tc cuprates and the k-(ET) compounds. However, the 
nodes exist at kz = ±7r and k^ ~ 0, therefore the SC or- 
der parameter is (ia;y-wave-like judging from the position 
of the nodes. Note that the b-c coordinate in the figure 
is square and is different from that of the original crystal 
for simplicity. 

Here we discuss the physical meanings of the calcu- 
lated results. The obtained Tn increases with increas- 
ing pressure up to 6 GPa. This behavior is totally dif- 
ferent from the result obtained by the mean field ap- 
proximation, which predicts a high Tn when the nesting 
condition is good and when DOS at the Fermi surface 
for U = {) (i.e., the inverse of the bandwidth) is large. 
This discrepancy comes from the fact that DOS at the 
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Fig. 3. ipi^yy^zy™T) on the Fermi surface at U=0.30 eV, 
T=14.4 K and the pressure of 15.2 GPa. An inset shows Fermi 
surfaces of U=0 (a dotted line) and U=0.30 eV (a solid line) at 
15.2 GPa. The unit of U is eV. 

Fermi level for finite U is reduced due to the singular to 
and /c-dependence of the self-energy 'E,k{uj) given by the 
FLEX approximation. The reduction in DOS is promi- 
nent as the Fermi surface becomes one-dimensional-like, 
because the singularity of Sk(t^) due to spin fluctuations 
becomes prominent in lower-dimensional systems. In 
other words, the increase in Tn is a natural consequence 
of the dimensional increase as stressed in the previous pa- 
per.^) As a result, the fine nesting does not always give 
higher Tn contrary to the present of the spin density wave 
in weakly correlated systems because strong fiuctuations 
in the lower-dimensional systems, or at lower pressures 
in this case, suppress growth on the long-range order. 
This is natural and plausible, and one of the merits of 
the FLEX approximation. 

Another important effect of the self-energy is that the 
nesting vector tends to be pinned to the commensurate 
value. The nesting vector at C/ = is not (tt, 0), however, 
it is (tt, 0) for moderate U owing to the k-dependence of 
Sfc(O) by the FLEX approximation. This effect naturally 
connects the magnetic structure in the strongly corre- 
lated (Mott) insulators with that in the strongly cor- 
related metals. The former is always commensurate, on 
the other hand, the latter can be incommensurate, but 
it becomes commensurate because of the strong electron 
correlations. It is also connected with the distortion of 
the Fermi surface at finite U; the electron correlations 
tend to make the Fermi surface nest better. 

With further increase in pressure, the nesting condi- 
tion deteriorates and the nesting vector becomes incom- 
mensurate. Then, Tn decreases and finally the antifer- 
romagnetic phase gets totally suppressed with increas- 
ing pressure. The SC phase that appeared in the high- 
pressure region is next to the antiferromagnetic phase 
and its SC order parameter has d-symmetry. The origin 
of the superconductivity is considered to be the AF fluc- 
tuation similar to those of (TMTSF)2X and k-(ET)2X. 
As a result, the obtained phase diagram. Fig. 2, as well 
as the concept of the superconductivity driven by the di- 
mensional crossover, will make sense and be reasonable. 

When one compares Fig. 2 with the experimental 
phase diagram, one notices that the pressure where the 
antiferromagnetism disappears is much larger than the 



4 J. Phys. Soc. Jpn. 



Letter 



Hiori Kino, Hiroshi Kontani, and Tsuyoshi Miyazaki 



experimental one. Here we enumerate the possible origins 
of this discrepancy: the error of the first-principles cal- 
culation, the approximation to calculate the self-energy 
and the model itself. 

The accuracy of the first-principles calculations for 
the title compound is investigated in ref. 10. Although 
it is well known that first-principles calculations based 
on the generalized gradient approximation cannot repro- 
duce the van der Waals interactions correctly, the in- 
vestigation shows that the error of calculated pressure is 
approximately 0.6 Gpa, which will not change the results 
shown in the present work. Similar results are reported 
for molecular solids and other charge-transfer molecular 
solids.15.16) 

Though it should be experimentally confirmed, we ex- 
pect that the first-principles calculation can determine 
the crystal and electronic structures of the title com- 
pound fairly accurately. Thus, one of the origins of the 
discrepancy on the transition pressure might be traced 
to the FLEX approximation. In a similar calculation 
of a nearly one-dimensional case at half-filling such as 
TMTSF salts, the SC phase is almost concealed with the 
AF phase in the FLEX approximation.^^) The particle- 
particle scattering channel or the vertex correction for 
susceptibility, which is not taken into account in the 
present approximation, may decrease the excessive AF 
region. 

We also comment on several problems concerning the 
efi'ective model. We fixed U in the obtained phase dia- 
gram for simplicity, however, such an assumption is vi- 
olated in real systems because the value of t{pl) con- 
siderably increases between P=0 and P=4 GPa. Al- 
though t(pl) is approximately constant for higher pres- 
sures, another problem emerges. The absolute value of 
t{p2) is considerably enhanced and becomes compatible 
with |t(pl)t under high pressures, which means that the 
assumption of the dimcr model, |t(pl)| ^ |t(p2)|, dete- 
riorates at a higher-pressure region. Therefore, i7eff at 
higher pressures is probably overestimated. As a result, 
[/eff, which is reasonably estimated to be 0.36 eV at 
ambient pressure from the value of t{pl), increases with 
pressure below 4 GPa. J7eff 0.6 eV at 4 GPa, and it de- 
creases for much higher pressures. This effect will make 
Tn and Tc smaller and reduce the critical pressure for 
the superconductivity. In addition, the dimer Hubbard 
model overestimates the spin fluctuation somehow, which 
will stabilize the AF phase too much. By solving these 
problems, we could obtain a theoretical phase diagram 
closer to the experimental one. 

Finally let us consider the reason why this system has 
a high Tc. The bandwidth increases almost linearly with 
increasing pressure as shown in Fig. 1. However DOS 
at the Fermi level does not decrease reciprocally as a 
function of pressure. There exists a saddle point near the 
r point that causes the peak (van-Hove singularity) in 
DOS. This peak shifts upward with increasing pressure. 



and the tail of the peak enlarges DOS near the chemical 
potential. As a result, DOS at the Fermi level at U=Q 
is almost constant, thus this system has a larger DOS 
at a chemical potential than the one expected from the 
bandwidth for P >4 GPa, which includes the pressure 
of experimental SC and theoretical SC. This finding is 
related to the high experimental T^. We note that in 
a quasi-one-dimensional system free from the van-Hove 
singularity near the Fermi level is very low (« IK), 
which is recognized by the FLEX approximation.^^) 

In summary, we studied the origin of the supercon- 
ductivity in /3'-(ET)2lCl2 based on the FLEX approxi- 
mation, using the tight binding fit to the electronic struc- 
ture in the first-principles calculation at ambient and 
under applied pressures. According to the present anal- 
ysis, a wide AF insulating phase is realized in the lower 
pressures region. Its magnetic ordering vector is approxi- 
mately (tt, 0). For the higher pressure region, next to the 
AF insulating phase, the dxy-wave superconductivity oc- 
curs owing to the AF fluctuations approximately (tt, 0). 
High Tc in /3'-(ET)2lCl2 is realized owing to the one- 
dimensional to two-dimensional crossover as a result of 
the crystal structure change under high pressures. Exper- 
imental checks for the magnetic ordering vector as well 
as the symmetry of the SC order parameter are highly 
demanded. 
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